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ABSTRACT 

In stellar evolution calculations, the local pressure scale height is often used 
to empirically constrain the amount of convective core overshoot. However, this 
method brings unsatisfactory results for low-mass stars (< 1.1 -1.2 Mq for Z= 
Zq ) which have very small cores or no convective core at all. Following Rox- 
burgh's integral constraint, we implemented an upper limit of overshoot within 
the conventional method of a parameterization in order to remove an overly large 
overshoot effect on low-mass stars. The erroneously large effect of core overshoot 
due to the failure of a parameterization can be effectively corrected by limiting 
the amount of overshoot to < 15 % of the core radius. 15 % of the core radius 
would be a proper limit of overshoot which can be implemented in a stellar evo- 
lution code for intermediate to low mass stars. The temperature structure of 
the overshoot region does not play a crucial role in stellar evolution since this 
transition region is very thin. 

Subject headings: convection — stars: evolution — stars: Hertzsprung-Russell 
diagram 



1. Introduction 

Understanding the physics of convective core overshoot is important in interpreting the 
color-magnitude diagrams (CMD's) and the luminosity functions of open star clusters, and 
generally of intermediate-age stellar populations. This is due to the facts that increasing the 
mixed core mass modifies the evolutionary tracks of stars (affecting the isochrone shape), 
and lengthens the evolutionary lifetimes of stars (affecting age determinations and luminosity 
functions near the main sequence turnoff). Thus, in addition to the classical problem of 



- 2 - 



establishing the chronology of star clusters in the Galaxy and nearby systems, core overshoot 
affects the spectral dating of young and intermediate age galaxies observed at large redshifts. 

Convective core overshoot is loosely understood here as the presence of material motions 
and/or mixing beyond the formal boundary for convection set by the classical Schwarzschild 
criterion (1906). An early study by Saslaw & Schwarzschild (1965), based on thermodynamic 
grounds (i.e. the edge of the convective core in massive stars is sharply defined in an entropy 
diagram), suggested that little overshoot takes place at the edge of convective cores. For this 
reason, it was believed that the presence of a gap in the CMD of open star clusters and it's 
magnitude could be used as indicators of the age and chemical composition of the cluster 
(Aizenman, Dcmarque & Miller 1966). However, Shaviv & Salpeter (1973) pointed out that 
if one takes into account the presence of hydrodynamic motions and turbulence, one might 
expect some non-negligible amount of overshoot. 

In fact, observational studies of the size of gaps near the turnoff in open star cluster 
CMD's suggest better agreement with theoretical isochrones which admit some amount of 
core overshoot (Maeder & MermiUiod 1981; Stothers & Chin 1991; Carraro et al. 1993; 
Daniel et al. 1994; Demarque, Sarajedini & Guo 1994; Kozhurina-Platais et al. 1997; 
Nordstrom, Andersen & Andersen 1997). Several computational schemes of various degrees 
of sophistication in treating the physics of overshoot have been developed for stellar evolution 
codes (Prather & Demarque 1974; Cogan 1975; Maeder 1975; Maeder & Mcynct 1988; 
Bertelli et al. 1990). As well, studies of detached eclipsing binaries in which the components 
have convective cores suggest some core overshoot of the order of 0.2 pressure scale height 
(Ribas, Jordi & Gimenez 2000). 

Stothers & Chin (1991) pointed out the high sensitivity of the amount of convective 
core overshoot to the adopted opacities. For example, increases in radiative opacities from 
the OPAL group (Iglesias & Rogers 1996) over the previous generation Los Alamos Opac- 
ity Library (Huebner et al. 1977), decrease the need for overshoot in comparison with 
observational data. In fact, the nature of the overshoot depends on the details of the local 
physics at the convective core edge. As pointed out by Zahn (1991), the local Peclet number, 
which characterizes the relative importance of radiative and turbulent diffusivity, determines 
whether in the mixed overshoot region, penetration takes place (i.e. the temperature gra- 
dient is the adiabatic gradient), or rather the local radiative transfer dominates the energy 
transport (overmixing) . In the latter case, the stable temperature gradient is unaffected by 
the mixing. 

Roxburgh (1989, 1992) also considered the physics of convective core overshoot from 
a different point of view. His integral constraint argument places an upper limit on the 
amount of core overshoot which can take place in the stellar interior. Zahn (1991) argued 
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that overshoot is practically adiabatic at the edge of the convective core, except in a thin 
transition layer, and therefore the Roxburgh's integral constraint defines the convective core 
size. 

More information can be gained from numerical simulations of convection. Because of 
the long relaxation times involved, it is not possible to perform 3D simulations of convective 
cores in full physical detail. Useful scaling information can however be derived from idealized 
3D simulations (Singh, Roxburgh & Chan 1995, 1998). 

The likely presence of differential rotation deep within stars brings further complexity 
to the overshoot problem, as one might expect a shear layer to develop at the edge of a 
convective core (Pinsonneault et al. 1991). Deupree (1998, 2000) considered the combined 
effects of rotation and convective overshoot in massive stars, which have large convective 
cores, with the help of 2D hydrodynamic calculations. Rotation likely plays a role in the 
case of less massive stars as well (in the range 1.5-2.0 Mq), with shear driven turbulence 
near the edge of the convective core inducing mixing from the helium enriched core into the 
envelope. 

Since a complete physical theory predicting the amount of overshoot for a star of a given 
mass and chemical composition, free of arbitrary parameters, is not available at this point, 
we approach the problem of core overshoot in the spirit of the work of Rosvick & VandenBerg 
(1998) in this paper. We consider and test parameterizations of convective overshoot, which 
are compatible with sound physical principles and with the growing data available on open 
cluster CMDs. These simple parameterizations have the merit of being readily implemented 
in a stellar evolution code. 



2. Test of Overshoot Treatments 

The local pressure scale height is often used to empirically constrain the amount of 
convective core overshoot in stellar evolution calculations. This conventional method assumes 
that the radial size of the overshoot region is proportional to the pressure scale height at the 
edge of a convective core, i.e. dos —a Hp. 

Observational studies of open clusters show that a moderate amount of core overshoot 
is essential to explain the observed gap near the turnoff in CMDs. In a study of NGC 2420, 
Demarquc ct al. (1994) found that a =0.23 is required for the best fit CMD. Also, Kozhurina- 
Platais et al. (1997) estimated a =0.20 and 0.25 for NGC 3680 and NGC 752 repectively 
by isochrone fitting (cf. Nordstrom et al. 1997 for NGC 3680). Overshoot effects appear to 
decrease due to the decreasing mass of stars near the turnoff in older clusters (Carraro et 
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al. 1994; Sarajedini et al. 1999). Although open clusters do not have many member stars, 
these studies of near solar metallicity clusters show that the conventional method constrains 
the amount of core overshoot relatively well for intermediate-mass stars. 

However, this pressure scale height method brings unsatisfactory results for low-mass 
stars (< 1.1- 1.2 Mofor Z= Zq). These stars have very small cores or no convective core 
at all, thus the local pressure scale height at the convective core edge would be very large 
implying a large amount of overshoot. Figure 1 illustrates the effects of the convective core 
and core overshoot on stellar evolution. With no overshoot, models of mass > 1.2 Mq show a 
MS hook which is caused by abrupt contraction of a convective core due to the flat hydrogen 
profile. When overshoot is included, evolutionary tracks show a redder MS hook due to 
the hydrogen supply from mixing in the overshoot region which lengthens the hydrogen 
core burning lifetime. However, the effect of overshoot on low-mass stars seems too large 
compared with observations. If the amount of overshoot is set to 0.2 Hp, low-mass stars with 
1.0 and 1.1 M0show a prominent overshoot effect in their evolutionary shape. However, 
these low-mass stars are not expected to have large overshoot regions since relatively old 
open clusters do not show an overshoot gap near the turnoff. Therefore, simple use of the 
a parameter cannot give consistent results for a wide range of masses and ages. 

Roxburgh's integral constraint (Roxburgh 1989; Zahn 1991) is another way to quantify 
the amount of core overshoot since it provides an upper limit to the extent of convective 
penetration in the following formula: 



where L^t is the total luminosity produced by nuclear reactions, L^afi is the radiative 
luminosity, and Vc is the radius of the effective core. Since the energy is transported by both 
radiation and convection inside a convective core, the integrand is positive until it reaches 
the Schwarzschild boundary where = L,.ad- The region of convective overshoot is located 
beyond this boundary where Lrad > Ltot up to the point r = Tc, which satisfies the constraint. 

However, Canuto (1997) pointed out in his theoretical investigation that Roxburgh's 
integral only gives an upper limit to the overshoot extent. Rosvick & VandenBerg (1998) 
used this constraint on NGC 6819 and found that a factor of 0.5 in Roxburgh's integral is 
required for the best fit CMD. Therefore, simply adopting Roxburgh's integral to quantify the 
amount of overshoot is not valid, and the fuzzy factor (analogous to the a parameter in the 
pressure scale height method) should be tested for clusters of various ages and metallicities. 

Faced with the current situation of the core overshoot treatments, we developed a con- 
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sistent approach which can be used regardless of stellar mass. Our simple approach of core 
overshoot is discussed in the next section. All stellar evolutionary tracks are constructed 
with Yale Rotating Evolution Code (YREC) with OPAL opacities (Iglesias & Rogers 1996). 
The general description of stellar evolution models can be found in Yi et al. (2001). 

2.1. Maximum amount of core overshoot (j3 limit) 

Roxburgh (1992) showed that his integral constraint gives an upper limit on the extent 
of convective penetration {dmax =P Rcore)- For very small convcctivc cores, the maximum 
amount of the extent of penetration is ~ 18 % of the core radius, independent of the details 
of energy generation and opacity. The maximum amount, P limit, varies from 0.18 to 0.4 
depending on the size of the convective core. 

Following Roxburgh's integral constraint, we implemented an upper limit of overshoot 
within the conventional pressure scale height method in order to remove the erroneously 
large overshoot effect on low-mass stars. In other words, we limited the core overshoot in 
a way that the amount of overshoot, dos (= a Hp), cannot be larger than a portion of the 
core radius, dmax (= P Rcore)- The merit of this modification is that the ad hoc size of 
overshoot regions for low- mass (< 1 Mq ) stars can be effectively corrected without affecting 
the evolution of higher mass stars. 

Since the integral constraint is obtained for convective penetration which maintains a 
nearly adiabatic temperature structure, we used the adiabatic temperature gradient in the 
overshoot region. The effect of the temperature structure will be discussed in the follow- 
ing section. Accepting the result of open cluster studies, we used a =0.2, and performed 
numerical experiments with various /3 limits for solar metallicity stars. 

The effect of the P hmit on stellar evolution is presented in Fig. 2. For low-mass 
stars, the presence of the overshoot effect is very sensitive to the choice of the upper limit of 
overshoot. It can be easily noticed that the overshoot effect on low-mass stars (mass < 1.1 
Mq ) becomes negligible for P < 0.15. However, if the upper limit is greater than 20 % of the 
core radius, the overshoot effect is still too large for low-mass stars. This result is consitent 
with Roxburgh's value of 0.18 for a nearly zero convective core radius. Thus, the erroneously 
large effect of core overshoot due to the failure of a parameterization for low-mass stars can 
be effectively corrected by limiting the amount of overshoot to < 15 % of the core radius. 

The choice of P does not significantly affect the evolution of intermediate-mass stars 
(1.5 - 2.0 Mq ) when /9 > 0.15 is used. These models are consistent with those of the same 
Q! without limiting the overshoot extent. However, if the amount of overshoot is restricted 
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to much less than 15 % of the core radius, the overall overshoot effect would be significantly 
reduced. Therefore /3 <S 0.15 for intermediate- mass stars gives effectively the same result 
with a smaller a of the conventional method, which is undesirable in our effort to modify the 
a paramerterizaion with the /? limit for a given a value. Combining the criteria for low-mass 
and intermediate-mass stars, we conclude that 15 % of the core radius (/3 =0.15) is a proper 
limit of overshoot for intermediate to low-mass stars. 

Due to the effects of overshoot on stellar evolution, such as shape of evolutionary locus 
near turnoff and time scale in hydrogen-burning stage, the ensuing isochrones are expected 
to show a different shape. We used 1.0 - 2.0 Mq evolutionary tracks with mass bin size 
0.1 in order to generate isochrones and test the effect of the p limit. We compared 
a— 0.2 isochrones with/without limiting the amount of overshoot (Fig. 3). As expected 
from the comparison of evolutionary tracks, (3 —0.2 isochrones are not much different from 
the isochrones without the P hmit. When overshoot is hmited to 10 % of the core radius, 
overall isochrones are dimmed due to the reduced size of the overshoot region, especially for 
younger ages. This result is similar to when smaller a is used. Thus, (3 =0.1 is too small to 
keep the amount of overshoot of intermediate-mass stars given by a =0.2. 

Focusing on the modification of the a parameter method, we increased /3 with increasing 
stellar mass in order to minimize the effect of limiting the extent of overshoot for stars of 
mass ^1.5 M© . Since the maximum extent of core overshoot depends on the size of a 
convective core, higher mass stars would have larger limits according to Roxburgh's integral 
constraint. Thus, we assumed a simple linear relation between stellar mass and the upper 
limit of overshoot {(3 = aM + b). The coefficients are determined for /3 = 0.1 and 0.4 for 1 
and 2 M© stars respectively. Fig. 4 demonstrates how the mass-/3 relation changes stellar 
evolution. Compared with (3—0.15 models, models with increasing (3 show slightly bluer and 
brighter isochrones for younger ages. Although increased P models are closer to the models 
without the /3 limit, the difference is not significant. 

2.2. Penetration and Overmixing 

The overshoot region generally represents a well mixed zone between convective and 
radiative zones. The temperature structure of the overshoot region depends on the physical 
conditions such as the ratio of radiative to kinetic timescales at the convective core boundary. 
According to Zahn (1991; also see Demarque et al. 1994), if the adiabatic cells penetrate into 
the radiative zone sufficiently far, then the temperature gradient of the overshooting zone 
would be adiabatic also (convective penetration). Or, if the adiabatic cells dissolve rapidly, 
then only mixing would occur with the remaining radiative gradient (overmixing). 
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However, it is not well known whether the temperature structure of the overshoot re- 
gion is adiabatic or radiative. In order to test the effects of the temperature structure on 
stellar evolution, we generated evolutionary tracks and isochrones with/without forcing the 
temperature gradient of the overshoot region to be adiabatic. 

Stellar evolutionary tracks of solar metallicity show almost no difference between over- 
mixing and convective penetration when the conventional method of a =0.1 - 0.3 is used. 
This is also true when the amount of overshoot is limited by any factor of the core radius 
(Fig. 5). We conclude that temperature structure of the overshoot region does not play a 
crucial role in stellar evolution, at least when the local pressure scale height method is used. 
This is due to the fact that this transition region is very thin. Note that the effect is slightly 
larger in low-mass stars. 

3. Conclusion 

This paper presents an approach to the treatment of convective core overshoot in stellar 
evolutionary models for stars with convective cores. By limiting the amount of convective 
core overshoot, the failure of the pressure scale height method can be effectively corrected. 
15 % of the core radius would be a proper limit of core overshoot for intermediate to low mass 
stars. We intend to use evolutionary tracks based on this approach to explore population 
synthesis and the spectral dating of high redshift galaxies. 

It is crucial to test stellar evolution by observation. Several observational approaches 
to test stellar models with convective overshoot are under way. The systematic study of 
open star clusters that exhibit hydrogen exhaustion phases provides important constraints 
that stellar models must satisfy. The WIYN Open Cluster Study (WOCS) project, pro- 
vides a comprehensive study of open clusters in the Galaxy selected to provide a range 
of chemical compositions and ages. In addition to photometric and spectroscopic observa- 
tions of individual stars, WOCS includes an astromctric proper motion study to establish 
cluster membership, a necessary condition. Because the number of stars in open clusters is 
small, establishing membership is essential for carrying out detailed tests of stellar evolution- 
ary models using a cluster CMD and luminosity function. This advantage is demonstrated 
in recent deep proper motion studies of NGC 188 (Dinescu et al. 1996) and NGC 3680 
(Kozhurina-Platais et al. 1997). An excellent overview of WOCS objectives and results 
obtained so far has been written by Mathieu (2000). 

Observations of the CMDs of some of the more populous star clusters in the Magel- 
lanic Clouds will provide more information on core overshoot in intermediate age stellar 
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populations. Although more distant than open clusters in the Galaxy, the Magellanic Cloud 
clusters have the advantage of testing the evolution of stars with a lower metallicity. 

Finally, stellar seismology, by analyzing oscillations that propagate deep in stellar in- 
teriors, promises to probe sensitively the convective cores of stars evolving off the main 
sequence in the intermediate mass range. Similarly, the mass of the exhausted helium core 
in subgiants and giants can be derived on the basis of physical models. Already observations 
of oscillations of nearby stars near the main sequence and in the subgiant and giant evolu- 
tionary phases have been done from the ground (Kjeldsen et al. 1995; Martic et al. 1999; 
Bedding et al. 2001) and in space (Buzasi et al. 2000). Preliminary analysis of some of these 
observations with the help of stellar models show the unmistakable signature (mixed p- and 
g-modes) of exhausted core size in the case of the more evolved stars (Christensen-Dalsgaard, 
Bedding & Kjeldsen 1995; Guenther & Demarque 1996). 
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Fig. 1. — The effect of convective core overshoot on stellar evolution for solar metaUicity 
stars of 1.5 to 1.0 Mq . Overshoot tracks have a redder MS hook due to the lengthened H 
core burning lifetime. The MS hook of low-mass stars in overshoot models show that the 
a parameter method fails due to the very large pressure scale height at the edge of the small 
convective core. 
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Fig. 2. — The effect of (3 limit on stellar evoluton for Z=0.018. The amount of core overshoot 
is limited to ^ factor of Rcore ■ Note that models with /? > 0.15 show nearly identical tracks 
for intermediate-mass stars (top panel). However, the overshoot effect is not present for low- 
mass stars when /? < 0.15 (bottom panel), which suggests that the erroneously large effect 
of overshoot on low-mass stars can be corrected by limiting the extent of core overshoot. 
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Fig. 3. — The effect of the (3 hmit on the isochrones of Z=0.018. /3=0.2 models are not 
much different from the models without limiting core overshoot, and repeat the failure of 
q; parameterization (top panel). When overshoot is limited to 10 % of the core radius, overall 
isochrones are dimmed, which is similar to the result of smaller a (bottom panel). 
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Fig. 4. — (3 limit as a function of the stellar mass. (3 is set up to be 0.4 for 2 Mq and 0.1 for 
1 Mq stars respectively. Compared with constant P models, models with increasing l3 show 
slightly brighter isochrones for younger ages. However, the difference is not siginificant. 
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Fig. 5. — The effect of the temperature structure of the overshoot region. The temperature 
structure of the overshoot region does not affect stellar evolution significantly since this 
transition region is very thin. 



